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FOREWORD 


This report Is a technical summary presenting the results of a 
study by the Electronic Systems Laboratory of the Electrical Engineering 
Department, Auburn University under the auspices of the Engineering 
Experiment Station. The report is submitted in partial fulfillment 
of the requirements of Contract NAS8-26752. 
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A STUDY OF POWER CONDITIONING 
AMD POWER DISTRIBUTION AMD COMPONENTS 

M. A, Honnell and H. M. Horton 
ABSTRACT 

The results of a comprehensive simulation and analysis performed 
on the operation of the regulator part of the Charger /Battery/ Regulator 
Module (CBRM) are detailed in this study. The CBKM is utilized as an 
integral component of the Skylab/Apollo Telescope Mount (ATM) electrical 
power system and contains a switching mode electronic regulator. Imp- 
lementing circuit analysis techniques* pertinent voltages and currents 
are calculated; these* in turn* are Incorporated into the regulator sys- 
tem study. Investigation of the turn-on and turn-off times associated 
with the switching circuitry is performed and an examination is made on 
these calculations* 

A computer program has been developed to provide a simulation of 
the switching regulator; this modeling is accomplished using the Fortran 
language and is processed on the IBM 360/50 computer. The inherent flex- 
ibility of the program allows changes in not only the input voltage and 
load resistance* but the important circuit component values as well. 

The simulation model computer program is utilized to generate 
graphs that will relate various CBRM parameters to one another. Through 
such graphs* the program will investigate the effects on output voltage* 
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output power » switching frequencies , and timing pulse-widths as input 
voltage and load resistance conditions are varied. 

The Sky lab /ATM has a complex electrical power system network; 
eighteen CBKM's are electrically connected in parallel and administered 
by a power sharing circuit. In order to analyze this parallel power- 
sharing arrangement, control techniques were employed and a feedback 
diagram was developed. Utilizing this feedback diagram, the system 
output voltage and individual regulator currents could now be formulated 
in tenaes of load resistances, distribution line resistances, and two 
intrinsic properties of the individual regulators: the output resistance 
and the no-load output voltage. 
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I. INTRODUCTION 


This report represents the results of a study made on the 
switching regulator of the Charger /Battery /Regulator Module (CBRM) 
which is a part of the Skylab /Apollo Telescope Mount (AIM) electrical 
power system. Eighteen such regulators , which supply common output 
buses, exist in the ATM and receive their energy from banks of solar 
cells. Basically, the regulator operates on the principle that a con- 
stant output voltage can be obtained at the output of a low-pass filter, 
if energy is correctly supplied to its input in pulses rather than con- 
tinuously. The majority of switching regulators come under one of two 
classes of timing control: 1) those having a constant on pulse-width 

and a variable switching repetition rate; and 2) those having a constant 
switching repetition rate and a variable on pulse-width. In the CBRM 
regulator both the on pulse-width and the switching repetition rate are 
variable. 

To provide the reader with an understanding of the CBRM's over- 
all operation, a simplified functional block diagram is presented and 
discussed in Chapter II. Also discussed are four basic possible con- 
figurations of the output stage of a switching regulator. Chapter III 
contains the development of mathematical expressions for the currents 
and voltages associated with the regulator and expressions for the 
switching times which produce these currents and voltages. 
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Using the information obtained in Chapter III, a computer pro- 
gram model of the CBBM regulator was developed and is presented in Chap- 
ter IV. Also included in that chapter is a discussion of the structure 
or format of the program and the first two pages of a typical output. 

The following are four features considered desirable in the computer 
program and consequently served as a guide during its construction. 

1) Furnish all voltages and currents pertaining to the output 
filter. 

2) Furnish output switching transistors' on time and dwell 
time, and furnish discharge time of the charging inductor. 

3) Provide a means of analyzing the load sharing properties 
of the regulator. 

4) Be versatile in format to allow changes in component values 
and input parameters. 

The problem of regulation and the reaction of the regulator to 
varying resistive loads are discussed in Chapter V. Also discussed are 
the nonlinearities of the regulator and the difficulties encountered 
in attempting to model such a system for purposes of analysis. 

Following the analysis of a single regulator. Chapter VI presents 
an investigation of the complete power system consisting of eighteen regu- 
lators in parallel and under the control of a power sharing circuit. To 
study the parallel arrangement, a feedback diagram is developed and the 
system output voltage and individual regulator currents are formulated 
in terms of load resistance, distribution line resistances and two in- 
trinsic properties of the individual regulators: the output resistance 
and the no-load output voltage. 
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Chapter VII contains a presentation of some of the results 
obtainable from the computer program model. The results are in the form 
of graphs displaying changes in one or more variables for a variety of 
input voltage and load resistance conditions. The conclusions given in 
this chapter summarize the work done, the results obtained, and suggest 
areas that could be investigated further. 



II, SYNOPTIC VIEW 


A discussion of the principles of operation of the regulator of 
the CBRM will be presented in this chapter and a detailed analysis of 
each part of the regulator will be presented in Chapter III. A simpli- 
fied block diagram of the regulator, showing the major blocks is shown 
in Figure 2-1* 

The principle function of the regulator is to take energy from 
one voltage supply whose voltage may be any value over a wide range (bank 
of solar cells) and use this energy to maintain an output voltage at a 
relatively fixed value. Since it is possible for the output voltage of 
the source supply to be less than the required regulated output, it is 
necessary to temporarily store this energy in the transition from input 
to output. The temporary storage element used is an inductor called the 
charging inductor located in the lower right-hand block of Figure 2-1 
and the energy is temporarily stored in the magnetic field of this induc- 
tor, The charging inductor is an integral part of an averaging filter 
which is also included in the same block of Figure 2-1. The energy from 
the inductor is transferred to large capacitors in the filter. The energy 
transferred from the supply source is equal to the energy supplied to the 
load connected to the output terminals plus a small amount lost in the 
switch and filter. 

When the switch in Figure 2-1 is closed, energy is transferred 
from the input supply to the magnetic field of the inductor and when the 
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switch is opened, this energy is transferred to the capacitors and load. 
Each time the switch closes and opens one increment of energy is trans- 
ferred from the supply source to the capacitors and load. The size of 
this increment, which must be dependent on the amount of energy being 
used in the load, is a function of the switch time and the input voltage. 
The on time of the switch is controlled by a one-shot multivibrator 
whose timing is controlled by other blocks in Figure 2-1. If the regu- 
lator is operated independently of other regulators, the switch time, in 
order to maintain constant output voltage, would be dependent only on 
the input voltage and the difference between the output voltage and re- 
quired voltage. However, the regulator in Figure 2-1 is designed to 
operate in parallel with several regulators to supply the total power 
to the load. So that all regulators furnish approximately the same 
amount of output power, their output currents are made the same. To 
accomplish this, the output current of each regulator is monitored and 
compared to a programmed current developed by a power sharing network. 
Each comparison produces a feedback signal which influences the timing 
of each switch to equalize the output currents. The output current mon- 
itor and power share blocks perform this function as will be explained 
later in Chapter VI. The output current monitor is also used to reduce 
the output voltage as a function of output current to give a character- 
istic as shown in Figure 2-2. The current monitor has a threshold value 
at 15.5 amperes that will cut the device down rapidly for overload pro- 
tection. 

In order to understand the basic operation associated with the 
incremental energy transfer from the input source to the output for the 
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Figure 2-2. The output volt age/output current characteristic. 
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purpose of voltage regulation, three basic circuits will be discussed. 
The first circuit has an output voltage that is controllable from zero 
up to the input voltage depending on the timing of the switch. The sec- 
ond circuit has an output voltage that is controllable from a greater 
value than the input voltage down to the input voltage, depending on the 
timing of the switch. The third circuit, has an output voltage that is 
controllable from a value greater than the input voltage down to zero. 

All three circuits will contain two switches, and S 2 » and 
will be analyzed under periodic operation of these switches. The period 
of the periodic operation is the time between consecutive closings of S^. 
The switch is assumed to be closed for seconds and opened for T 2 
seconds, thus giving a period of plus T 2 seconds. The second switch 
is synchronized with and is opened when is closed and closed 
when is opened. For convenience, a term D has been defined as the 
duty cycle of S-^ and is the fraction of the period during which is 
closed and is given by D * + T 2^* 

A fourth circuit to be discussed is a modification of the third 
example and is the actual circuit used in the regulator of the CBRM. 

Four Basic Configurations of Output Stages 
1. Buck Circuit 

The circuit in Figure 2-3 will have an output voltage whose 
average value will be between zero and V g depending on the duty cycle of 
switch S^. The voltage v(t) in Figure 2-3 is shown in Figure 2-4. A 
linear analysis of the circuit to the right of switch reveals that 
the output voltage will have the same average value as v(t). This 
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Figure 2-5. The range of possible values for the output voltage Vq 
a function of the duty cycle. 


as 
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analysis is based on representing v(t) as a Fourier series and realizing 
that a constant v(t) will produce a constant Vg(t) of the same value. 

The harmonics of v(t) will also appear in Vg(t) ; however » they will be 
attenuated by the low-pass LC filter. The average value of v(t) , and 
vg(t) is: 


<v(t)> - <v ( t )> - = DV S (2-1) 

X + 2 

and is plotted in Figure 2-5 as a function of D* 

II. Boost Circuit 

The circuit in Figure 2-6 is commonly known as a Flyback circuit 
and its output voltage will be between a value greater than V g and V g 
depending on the duty cycle of In order to solve for the output 

voltage during steady-state* steady-state conditions are assumed. It is 
also assumed that during the time that is closed* the voltage Vg(t) 
will not change an appreciable amount. Under these conditions the vol- 
tage v(t) in Figure 2-6 will have the waveform shown in Figure 2-7. 

Since the current in L is continuous and periodic, the average value of 
the voltage across L over one period will be 



(t)dt 



|ti(T) - i(0)J. 


( 2 - 2 ) 


Under steady-state conditions the value of the current through the induc- 
tor at the end of a period, T, is equal to the value of current through 
the inductor at the beginning of the period. Thus , by Equation 2-2 the 
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average value of the voltage across the Inductor is zero. Solving for 
the average value of v(t) , the voltage across S 2 : 


v(t) » V s - v L (t) 

(2-3) 

<v(t)> *■ <v s > - <v L (t)> 

(2-4) 

<v(t» <V g > 

(2-5) 

The average value of v(t) in terms of 
2-7: 

Vq can be calculated using Figure 

v(t) - ft V ° - DV 0 

1 + 2 

(2-6) 

Combining the two equations for v(t) 

and solving for VqJ 

V = — V 
0 D s * 

(2-7) 


Vq is plotted as a function of D in Figure 2-8. 

III. Buck-Boost Circuit 

The circuit in Figure 2-9 will have an output voltage whose 
average value will be between a value greater than the input voltage and 
zero, depending on the duty cycle of switch As was done for the 

Boost configuration, an expression for Vq as a function of D can be de- 
rived for steady-state conditions with a continuous and periodic current 
through L. Again the average value of v(t) , the voltage across the indue- 





16 


tor is zero. From Figure 2-10 which shows v(t) as a function of time, 


<v(t)> 


T l v 0 + T 2 v s 
T l + t 2 


0 . 


Solving for Vq from this equation yields: 


DVq + (1 - D)V S - 0 



(2-8) 


(2-9) 

(2-10) 


The possible range of Vq is shown in Figure 2-11 and as indicated by 
Equation 2-10 will have negative values. 

IV. Modified Buck-Boost Circuit 

The output stage used in the GBRM regulator is a modified version 
of the Buck-Boost configuration just discussed. As shown in Figure 2-12, 
a segment of time (known as dwell time) is introduced in which both 
switches, and S2» will be opened. A drawing of this configuration is 
given in Figure 2-13 and the waveform for v(t) appears in Figure 2-14. 

The average value of v(t) will be zero as discussed for the two 
previous configurations and with the aid of Figure 2-14 can be expressed 
as : 


/ V ( t )S * li v + h — 

\ ' T 1 + t 2 + t 3 ° T 1 + t 2 + T ' 


V g - 0 


( 2 - 11 ) 


DV q + 


T 1 + T 2 + T 3 - T i - t 3 

T 1 + t 2 + t 3 


0 


( 2 - 12 ) 
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\t- ONE PERIOD H 


Figure 2-12. A graph of the switch position relationship between SI 
and S2 of the Modified Buck-Boost output stage. 
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Figure 2-14. The waveform of v(t) , the voltage across L of Figure 2-13. 
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DV Q + (1 - D - D d )V s - 0 


02-13) 


_ 1 - D - p d 


02-14) 


where D d 


T 1 + T 2 + t 3 


(2-15) 


and Is the fraction of tine when both switches are opened, A plot of 
Vq as a function of D for a constant is shown in Figure 2-15, 

In the actual circuit of the CBRM regulator, large amounts of 
current must be switched by S. and S- and carried by L, To facilitate 

X * 

handling this current, four branches each composed of an S^, and L 
are paralleled and each takes a fourth of the total current. Switch 
is a silicon power diode and is a silicon power transistor. This 
arrangement is shown in Figure 2-16. For purpose of analysis, the four 
inductors in parallel can be considered as one with an equivalent value. 




transistor switches , four 




III. MATHEMATICAL DERIVATIONS 


Derivation of the mathematical equations governing the operation 
of the regulator will be developed based on the piece-wise linear behav- 
ior of the system. To understand the piecewise-linear behavior, let it 
be assumed that the regulator is in a normal mode of operation with an 
output current being delivered to a load. The output voltage is continu- 
ously compared with an internal reference voltage and if low-voltage 
condition is indicated, the one-shot multivibrator is triggered. The on 
time of the multivibrator is a function of the input voltage and the out- 
put current. The output of the multivibrator is used to turn on transis- 
tor switches which allow the charging inductor to receive energy from the 
input source. The amount of energy stored in the inductor is a function 
of the on time of the multivibrator and the value of the input voltage 
supply. The higher the input voltage the shorter the time required 
to store a given amount of energy. Also, since the energy delivered to. 
the load must balance with the energy transferred through the inductor, 
it is necessary that the on time of the multivibrator increase as a 
function of the load current. Therefore, the on time of the multivibra- 
tor is a decreasing function of the input voltage and an increasing 
function of the output current. When the multivibrator turns off the 
transistor switches, the energy in the charging inductor is transferred 
into the output filter as explained in Chapter II. When the current in 
the inductor reaches zero, the energy dump from the inductor to filter 
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is complete and control of the regulator is returned to the comparator. 
The multivibrator is again automatically triggered as soon as the com- 
parator determines that a low-voltage condition exists. This sequence 
of events constitutes one cycle of operation of the switching regulator. 

During steady-state operation, there is a time interval between 
the end of the energy dump and the condition necessary for the retrig- 
gering of the multivibrator. This is due to the building up of the out- 
put voltage by the energy just dumped by the charging inductor, and then 
waiting for this voltage to decrease. This wait period is usually called 
dwell, or off, time. 

One cycle of operation consists of three distinct time intervals; 
the inductor charging time, the inductor dumping time, and the dwell 
time. The operation of each part of the regulator will be analyzed for 
the three time intervals. The output filter is analyzed first followed 
by an analysis of the charging inductor. The results of these analyses 
will then be used to determine the necessary timing equations for satis- 
factory operation of the regulator. 

A. Output Filter Currents 

The output filter is shown in Figure 3-1 and is a fourth-order 
network when the commutating diode is forward biased (dumping time in- 
terval) and is a third-order network when the diode is reverse biased 
(charging and dwell time). The mathematical analysis is divided into 
two parts depending on the state of the diode. The mode of operation 
during the time interval when the diode is forward biased will be re- 
ferred to as mode A and the mode of operation during the time interval 


when the diode is reverse biased will be referred to as mode B. 
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The equations for the filter currents will now be derived for 
mode A operation. A matrix equation for the loop currents, in Laplace 
transform terms, is given on page 27 as Equation 3-1. In the matrix 
equation and the following analysis (for mode A), the diode is modeled 
as a small voltage source, V D , in series with a small resistance. The 
form of I A i(S), I^tS), and ^(S) is! 


! a (S) 


-v T 


1 

- + 


+ a^S + 


*D + «11 + *L2 + »L S V* + b 3 s3 + b 2 s2 + b l S + b 0 


(3-2) 


The derivation of this equation and expressions for the coefficients are 
given in Appendix A. 

The last term of Equation 3-2 has a denominator which can be 
expanded into two complex conjugate pairs of roots and can be written 
in the form of : 


last term 


R A1 / e Al r A 1 /" e Al R A2 / 0 A2 r A2 / -0 A2 

* * ! "■'— ■■■ ■ + ■■ + 

S+0t Al ” j W A1 S+ct *t + ^*i s+a a ^ “J w ^ S+a^+jw, 


(3-3) 


A1 J A1 


A2 A2 


A2 J A2 


Thus, the time domain expression for the currents will have the form of: 


i A (t) - 2R A1 expt-a^tj cos(w A1 t + 6^) 

+ 2R a2 expt-a^t] cos^^t + 9 A2 ) 

- V ( *D + *L1 + hi + V • (3-4) 

Each parameter in the current equations is defined and discussed in 
Appendix B. 



( ia + su + *li + “ci + V 


(- — - Rr, ) 

SCI C1 


SH - V 


+ rrr + su + R cl + Rc 2 + r l2 > (- ^ - R^) 


SCI SC2 


- w 

SC2 ^ 


<Ia + + *«> 




Ia 2< s > 


I A3 (S) 


LI • 1 L1 (0) - v a (0)/S - V D /S 


L2 • i l2 (0) + v cl (0)/S - v c2 (0)/S 


'hJ 


(3-1) 


V C2 (0>/S 


where 1 ^ 1 (0) 

i L 2<°> 

v c lC0) 

v c2 (0) 


is the initial current through 11, 
is the initial current through L2, 
is the initial voltage across Cl, 
is the initial voltage across C2, and 


Vjj is the forward voltage drop across the diode. 
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When the diode is reverse biased, it is in mode B operation and 
appears as shown in Figure 3-2. A matrix equation for the loop currents 
in Laplace transform terms is: 


^SCl + SC2 + SL2 +R C1 + ®C2 + ®L2^ SC2 " *C2^ 


- 

I B2^ S ^ 

SC2 * ^SC2 + R C2 + ®L^ 


Ib3<S) 


L2 * i u (0) + v cl (0)/S - v^COJ/S 


(3-5) 


V C2 <0)/S 


where i^(0) is the initial current through L2 
v cl <0) is the initial voltage across Cl, 
v C2 (0) is the initial voltage across C2. 
The form for Ig 2 (S) and I B3 (S) then will be: 


* 

and 


I,CS> 


CjS* 4* c^S + Cq 
d 3 S 3 + d 2 S 2 + djS + d Q 


(3-6) 


Equation 3-6 has a denominator Which can be expanded into a 
complex conjugate pair and one real root and can be written in the 


form of: 
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B. Model of Charging Inductor 

If the charging inductor, LI, is ideal and the voltage across it 
is a constant, then the current through it is a linear function of time 
and is as follows: 



(3-9) 

(3-10) 


where is the voltage across the inductor. 

However, for large currents, the inductor's value will be a function of 
the current through it. Under this condition, the current through the 
charging inductor will still be a function of time but not a linear one, 
and can be expressed as: 



(3-11) 
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Since the charging inductor's value does vary with current, it 
is necessary to model the inductor so that this variation is accounted 
for. Each charging inductor in the CBRM power supply has a value of in- 
ductance of approximately 200 pH when operated well below saturation. 

The voltage across LI can be expressed as: 



( 3 - 12 ) 


and holds for any set of time and current values when the., proper value of 
LI is used. If a very small change is taken in time, Ll(i) will be a 
near constant. Thus: 


V L = U(i, ^ 


(3-13) 


The smaller the At, the greater the accuracy of this equation. Solving 
for the change in current: 


Ai(t) 


VrAt 

11 ( 1 ) 


(3-14) 


This equation will give the increase in current through the inductor for 
a small increase in time when the inductor's value is known in the range 
of operation. 

Equation 3-14 states that for fixed increments of time, At, the 
increments of current, Ai, get larger as Ll(i) decreases. Figure 3-3, 
drawn from measurements on the actual CBRM regulator, shows the results 
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of such action. The input voltage and load current were 60.3 volts and 
14.6 amperes, respectively. 

It can be seen by Equation 3-14 that the inductance for any value 
of current is equal to the applied voltage divided by the slope of the 
current curve at the operating point. 


U(i) 



Ai(t)/At 


(3-15) 


From measurements on the regulator system, it was calculated that the 
value of the inductance at two amperes is approximately 200 pH and at 
twenty amperes the value is approximately 100 pH. 

One method of determining the current through the inductor is to 
use Equation 3-14 and represent Ll(i) empirically as: 

Ll(i) * 200 pH - Kl(i)^ (3-16) 


This equation was chosen because it generates values of LI close to ac- 
tual values and holds for all i within the range of operation. Kl and 
K2 are selected constants used to make the equation approximate the ac- 
tual inductor. Values used for Kl and K2 that gave good results were 
0.70(10)“® and 1,08, respectively. 

An iterative process will now be developed that will relate a 
value of current through the inductor, LI, to a given fraction of the 
total time for which a constant voltage is applied across the inductor. 
Let the total time be T and the applied voltage be V L . The time T is 
divided into small segments At. Beginning with an initial current of 
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zero, Ll(i) is calculated by Equation 3-16 and used in Equation 3-14 to 
determine the amount of increase in current, 61, due to an increase in 
time At. How using the new value of current, another value of Ll(i) is 
calculated and used to calculate the latest increase in current due to 
another increase in time. This procedure is continued until the sum of 
the At 's equals T, each time using the sum of the Ai's to determine 
Ll(i). This process yields very satisfactory results if At is T/20 or 
smaller. 

C. Dwell Time 

Whenever the voltage at the input of the output filter is suf- 
ficiently high, the one-shot multivibrator should stay off. This off 
time is called dwell time. Looking at the diagram of the comparator 
circuit in Figure 3-4, which turns the one-shot on, the voltage equation 
for its difference-mode input can be written: 


15 R3 - 7.5K V„, 

V S : — _ y 

a 7.5K+R3 CURHT 


(3-17) 


where is the voltage across the filter capacitor Cl, and is 

the voltage at the inverting input of IC2 produced by the output current 
monitor. The transition of from a negative value to a positive value 
will cause the comparator* s output to go from negative to positive there- 
by turning the one-shot on. 

In determining the dwell time, the question to be answered is 
,f How long does it take for the first capacitor, Cl, in the output filter 
to discharge to a value that would require the one-shot to come on again?" 
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has a waveform similar to that shown in Figure 3-5. The exact ampli- 
tudes and times are functions of the input voltage and the load current. 

Let the value of V_, at the end of the dump time, which will be 

Cl 

discussed later, be called V and let the value of V that makes V a 

Cla Cl a 

go from a negative value to a positive value (V ■ 0 + ) be called V 

a C1B 

Then the decrease in V . during the dwell time will be: 

Cl 


AV, 


Cl 


Cla 


- V 


CIS 


(3-18) 


To develop an expression for the dwell time, another equation for AV^ 
can be used if the current being discharged by Cl is considered a con- 
stant. Due to the large value of capacitance and short dwell time, this 
is a good approximation. Therefore representing dwell time by T : 



(3-19) 


AV 


Cl 


X B2 T D 


Cl 


(3-20) 


Combining Equations 3-18 and 3-20 and solving for V 

which V equals 0 + yields: 
a 


Clg* 


the voltage for 


IT 
B2 D 


V = V 
C16 Cla 


a 


(3-21) 




Figure 3-5. A typical waveform of V 


the voltage across Cl. 
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Using this value o£ Vq in Equation 3-17. and setting V a ■ 0 and solving 
for the dwell time, T D , yields: 


T 


D 


V CUENT (7,5K + ^3) - 15 R3 


+ V Clct 



( 3 - 22 ) 


Thus, when the voltage at the inverting terminal, the vol- 

tage across Cl at the end of the dump time, and Ig2* current in loop 
two of the output filter, are known, Tq can be calculated, can 

be determined through an analysis of the SC circuit connected to the 
inverting input terminal. 

D. Charging lime 

The charging time is the length of time that the one-shot multi- 
vibrator must stay on. The one-shot shown in Figure 3-6 is normally off 
with -13 volts at the output. Under this condition there is approximately 
-0.6 volt at the inverting input of IC1 due to CR8 being forward biased 
and there is approximately -1.7 volts at the noninverting input due to 
the R18, B23 voltage divider network (CR10 is forward biased and CR6 is 
reverse biased) . The one-shot remains in this stable condition until the 
current in the charging inductor is zero and the output of the comparator 
goes positive. When this happens, the noninverting input of 1C1 goes 
positive with respect to the inverting input, the one-shot is turned on, 
and the output of 1C1 goes to 15 volts positive. There will then be 
current through the charging inductor and when it is detected, the output 
of IC2 will be drawn to ground potential prohibiting the comparator from 
having control of the one-shot until the inductor has been charged and 
completely discharged into the filter. With the comparator's output 
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grounded, there is a positive voltage on the noninverting input of 1C1 
from the current monitor which is proportional to the output current. 
Since CR12 is now reverse biased, C6 will begin to charge . The time 
interval it takes for the voltage at the inverting input to just exceed 
that at the noninverting input is the on- tine for the switching trans- 
istors and charging time for LI, and will be referred to as T^. At the 
end of this time, the one-shot will turn off and wait for the next turn- 
on command. 

An expression for in equation form can be developed as fol- 
lows. First, an equivalent circuit for the input to the 709 operational 
amplifier, I Cl, and the charging path for C6 is modeled in Figure 3-7 
in which the input resistance to IC1 has been named R21. Solving for 
the time required to charge C6 to the value that makes the voltage 
across R21 equals zero will yield: 


„ (R12 + R20HB21 + R22) C6 - 
C R12 + R20 + R21 + R22 

(3-23) 


v in - \ 


where (R3 + R17 [{ E18) |j R23 has been named R22, V^ n is the input voltage 
and 


R23 (0.9 I B2 ) 

R23 + R3 + R17 jj R18 


(3-24) 


" 0.2 1 B2 


(3-25) 



RI2 + R20 



OUTPUT CURRENT 
MONITOR VOLTAGE 


Figure 3-7. An equivalent circuit of the one-shot multivibrator circuit shown in Figure 3-6 when 
it is in the on state. R21 represents the internal resistance of the operational 


5Tht>; 
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The voltage, 0.9 1 q 2> i- n Equation 3-24 is the voltage appearing at the 
output- current-monitor voltage terminal shown in Figures 3-6 and 3-7 
and is developed as follows. The current in the second loop of the out- 
put filter is monitored by a transductor which controls an oscillator. 
The output of the oscillator is rectified and filtered and produces 
0.1 volt per anpere of current. This signal is simplified by a factor of 
nine and then appears at the output -current-monitor voltage terminal. 

Note that if the input resistance, R21, is infinite, the equa- 
tion for T c reduces to 


T 


C 


(R12 + R20) 06 In 



0.6 





(3-26) 


But such is not the case. Typical input resistance for a 709 amplifier 
is 400 Kil and can be as low as 130 K0. This large but noninfinite value 
must be considered as an important factor because R12 is in the 600 Ki3 
to 900 KS1 range* 

E. Energy Dump Time 

The last interval of time to be determined is that needed for 
the charging inductor to dump. This energy dump time is called Tg and 
will be found with the aid of Figure 3-8 which shows only the input 
branch of the output filter. Due to the relative values of the filter 
components, most of the current being discharged by LI is conducted by 
the input branch. As shown, this branch contains a resistor, R^, to 
represent the small effective resistance and an inductor, , to 
represent the effective inductance of Cl. 
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+ V IN 



-V| N 


Figure 3-8, The input circuitry of the output filter stage; 

used to determine the dump time for inductor Ll 
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It was developed in section B covering the charging inductor, and 
indicated by Equation 3-14, that the amount of change in current through 
the inductor for a constant applied voltage, is a function of the length 
of application time. The voltage across the inductor in the present 
circumstances will be -(V^ + 0.6) volts and will not be constant through- 
out the dumping interval due to a small change in v a ; however, this does 
not affect the validity of the equation. Therefore, the change in current 
for a small change in time will be: 


Ai(t) = 


(V a + 0. 6)At 
Ll(i) 


(3-27) 


Solving for At: 



(3-28) 


A suitable Ai(t) can be chosen by dividing the initial value of the dump- 
ing current, which was the final value of the charging current and will 
be known, into small segments. 

In a manner similar to that used in determining the charging 
current, the dump time can be found by using an iterative process. After 
determining a Ai(t) , Equation 3-16 can be used to calculate Ll(i) for 
the peak value of the charging current. Using this Ll(i), a At can be 
calculated by Equation 3-28. An expression for during this time 
will be developed later. Adding Ai(t) , a negative number, to the init- 
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ial current will give the remaining current which is used to calculate 
the next At. Each time through the iterative process, the remaining 
current through the inductor is determined and used to calculate a new 
value for Ll(i) which in turn is used to calculate a new At. When the 
current has been incremented down to zero the sum of the At's will be 
the total dump time, Tg. 

The voltage, , across the branch containing Cl can now be 
calculated for the interval Tg. This voltage will be made up of a large 
constant voltage plus a small ripple voltage. Although its name does 
not explicitly indicate it, will be a function of time and as shown 
in Figure 3-9 can be expressed as: 


V C1 ■ »cl< c > + “ell 4 Ut> - 1*2 00] + V Lcl (3-29) 


where: VQ^(t) is the voltage across Cl only, 

i(t) is the current dumped by each charging inductor, and 
iA 2 ^t) is the current in the second loop of the filter* 
Since most of the current being dumped goes only into the branch contain- 
ing Cl and since the voltage across this branch as well as the output 
voltage is a near constant, *£ 2 ^) changes very little and will be con- 
sidered a constant named 

The total current being dumped is four times i(t) and equal to: 

4i(t) = 4 (mt + I D ) - Mt + 4 I 
f P 


(3-30) 
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where m is the slope of i(t) in amperes per second and Xp is the 
peak value of i(t) at t ** 0 + . 

Considering the slope a constant, v^(t) can be written as: 

[Mt + 4 I - I^jdt + Vq (0 + ) (3-31) 

where v { ,^(0 + ) is the initial voltage across Cl. 

Evaluating the integral yields: 

v ci (t> ■ -ff[f t2 + w h - W*] + «- 32 > 


'Cl 


(t) 


-JLf 

01 J 

Jr\ 


While M is not actually a constant it is near so and an expression for 
it can be found in the interative process developed to calculate Tg. 

In that process, a Ai(t) and At were found as the charging inductor 
dumped. Since four such inductors are present: 


At 


(3-33) 


Thus in the process of calculating the dump time Tg, can be calcu- 
lated. Each time a At is determined, it is used with the Ai(t) to cornr- 
pute by the following; 
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_l[ 2 iiitl t 2 + (4 Ip - tA2 )tj + Vcl (0 + ) 


+ [4 i(t) - I w ] + V Lcl . 


(3-34) 


Before the above equation can be applied , Vj^ ““St he known. Consider 
Figure 3-10 which approximates the current flowing through one transistor 
switch, the current being dumped into Cl from one charging inductor, and 
the current in one charging inductor. The time interval r 0 £f is the 
turn-off time of the transistor switches and has been exaggerated in 
the drawing. For the transistor switches in the regulator, will 

be less than or near, one microsecond depending on the base drive and 
the amount of current being switched. 

As long as the switches are in the process of closing, the 
total voltage across is : 


(3-35) 


if a constant slope for the current is assumed. Curing Tg, shown in 
Figure 3-10, the total voltage across Lq is equal to: 


V L C1 - L C1 


(3-36) 
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V L a “ 4Lcl ^ <3 - 37) 

and will have a negative value. In the computer program. Equation 3-35 
will be used to calculate Vj^ during the first microsecond after the 
charging time, Tq, and Equation 3-37 will be used thereafter. 

Using the equations developed in this chapter, it is now possi- 
ble to construct a computer program model of the GBRM, The output volt- 
age will be as shown previously in Figure 2-2 on page 7, and can be 
expressed as: 


YquT " 30*4 ~ 0.213 Iflyy 


(3-38) 


where 30.4 is the Thevenin voltage (no-load voltage), 0.213 is the 
Thevenin resistance (effective dc output resistance) , and Iq^ is the 


output current. 



IV. COMPUTER PROGRAM MODEL 

This chapter contains the listing of a computer program model 
of the CBEM regulator which was developed using the mathematical analy- 
sis done in Chapter III. For clarity, the program contains comment 
statements which serve as a guide to the reader. The discussion in the 
following section is given so that the structure of the program will be 
understood. 

A. Program Construction 

In the computer program, the switching times are assigned new 
descriptive names and the text of this chapter, as well as those follow- 
ing, will conform to the assignment which is as follows: 

Charging time (T c ) : TCHARG 

Dwell time <T D ) : TDWELL 

Energy dump time (Tg) : TDUMP 

Figure 4-1 shows a timing chart for the regulator and starting 
with a given set of initial conditions, the computer program must be 
able to simulate the cyclic process shown in this timing chart. Besides 
the load resistance R L , some other component values are left as vari- 
ables at the beginning of the program to allow changes to be made. 

These components are listed and identified at the beginning of the 
program. 
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Figure 4-1. Timing chart showing sequence of the switching times. 


Foe a given set of component values, the roots to the current 
equations (3-2 and 3-6) of the output filter are fixed. Since the roots 
are needed to determine the denominators and numerators of the partial- 
fraction-expansion forms of these current equations, they are calculated 
at the beginning. To do this, an IBM scientific subroutine called 
POLRT, for Polynomial RooTs, is used. All it requires are the coef- 
ficients of the denominators of Equations 3-2 and 3-6. Next in the 
program, the denominators of the partial fraction expansion are calcu- 
lated and as indicated in Equation 3-3 and 3-7 they can be complex. 

The above operations occupy blocks one through five in the pro- 
gram flowchart shown in Figure 4-2. It is intended that this chart 
show only the major functions to be performed and without details. 

Blocks six through nineteen represent the cyclic action of the 
regulator. The process begins with the output filter in mode B; the 
transistor switches are off and the charging inductor contains no 
current. The length of time that the transistor switches remain off, 
TDWELL, must be calculated. Once TEWELL is known, CURNT(2) can be cal- 
culated at the end of that time and this current can be used to calcu- 
late TCHABG. During TCHARG the transistor switches will be on. As 
shown in the timing chart, at the end of TCHARG the switches will turn 
off and the charging inductors will begin to dump, placing the filter 
in mode A. Thus all filter currents and voltages must be calculated 
and used as new initial conditions for equations relating to mode A. 

Also CHARGI, the current through the charging inductor must be calculated. 
The filter will stay in mode A as long as the inductor is dumping. At 
the end of the dump time, TDUMP, the latest values of filter currents 




Figure k-2. A simplified flowchart of the computer program. 
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and voltages are calculated and used as the Initial conditions going 
into mode B. The cycle has been completed. 

B, Program Listing 

A listing of the computer program begins on page 56. First 
there is a glossary of terms to assist in following the program. This 
is followed by a list of components and initial^ conditions. Following 
the program listing, are the first two pages of a typical printout. 

The computer source program is written in IBM System/360 
FORTRAN IV which is compatible with USA FORTRAN. As stated previously, 
one IBM scientific subroutine, FOLRT, is used. In addition, several 
built-in, FORTRAN -supplied mathematical function subprograms are used, 
such as Absolute value (ABS), Exponential (EXP), Natural logarithm 
(ALOG), and Sine (SIN). 

The computer used was an IBM System/360 Model 50. To produce 
the plots given in Chapter VII, arrays to store data points were placed 


within the program and their contents plotted on a Calcomp Plotter. 
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SWITCHING REGULATOR COMPUTER MODEL **** 


GLOSSARY OF TERMS USED IN THIS PROGRAM 

ALPHA4 o ALPHAS- THE DAMPING FACTORS® EQUAL TO MINUS THE 

REAL PART OF THE OUTPUT FILTER'S ROOTS. 
CHARGI- CURRENT IN THE CHARGING I NDUCTOR * LI. 

CKA,CXB- COEFFICIENTS OF PARTIAL FRACTION EXPANSION OF 
CHARACTERISTIC EQUATIONS. 

CNUMA* CNUMB- NUMERATORS OF CHARACTERISTIC EQUATIONS* 

USED IN FINOING CKA AND CKB. 

CURNT2- MAGNITUDE OF THE CURRENT THROUGH L2. 

C1INTV- INITIAL VOLTAGE ON Cl. 

C2INTV- INITIAL VOLTAGE ON C2. 

C3INTV- INITIAL VOLTAGE ON C3. 

CYCLE- RUNNING COUNT OF THE CYCLES THE REGULATOR HAS 
GONE THROUGH. 

DNOMA® 0NOMB- DENOMINATORS OF CHARACTERISTIC EQUATIONS* 
USED IN FINOING CKA AND CKB. 

Kl* K2- CONSTANTS USED IN MODELING SATURATION 
CHARACTERISTIC OF LI. 

L1INT1 - INITIAL CURRENT THROUGH LI. 

L2INTI- INITIAL CURRENT THROUGH L2. 

LIS- ACTUAL INDUCTANCE VALUE OF Li® IT WILL BE A 
FUNCTION OF THE CURRENT THROUGH IT. 

N- SUBSCRIPT FOR CALCULATED OATA POINTS* BEGINS AT 1. 
MAXC- MAXIMUM NUMBER OF CYCLES PROGRAMMER DESIRES THE 
REGULATOR TO GO THROUGH. 

OMEGAA* OMEGAS- THE OAMPEO NATURAL FREQUENCIES OF THE 

OUTPUT FILTER. EQUAL TO THE POSITIVE 
IMAGINARY PART OF THE OUTPUT FILTER'S 
ROOTS. 

RA* RB- MAGNITUDE OF CURRENTS IN OUTPUT FILTER. 

ROOTA, ROOTS* ROOTS OF CHARACTERISTIC EQUATIONS OF 
OUTPUT FRTER. 

T— TI ME . 

TCHARG- TIME FOR INPUT INDUCTOR* LI * TO CHARGE. 

THETAA 9 THETAB- PHASE ANGLE OF CURRENTS IN OUTPUT 

FILTER. 

TTOTAL- ACCUMULATED TIME* EQUALS SUM OF ALL TDWELLS* 
TCHARGS* AND TDUMPS. 

TDUMP- TIME FOR INPUT INDUCTOR* LI® TO DUMP. 

TDWELL- REGULATOR DWELL TIME. DURING THIS INTERVAL, THE 
TRANSISTOR SWITCHES ARE OFF. 

VC1- VOLTAGE ACROSS THE INPUT CAPACITOR* Cl* OF THE 
OUTPUT FILTER. 
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VCURNT- THE VOLTAGE AT THE NEGATIVE INPUT OF IC2-A16, 
A COMPATATOR CIRCUIT WHICH DETERMINES TDWELL. 
VIN- REGULATOR INPUT VOLTAGE. 

VOUT- REGULATOR OUTPUT VOLTAGE. 

DIMENSION COFI4I,XCOF<5I»ROOTII4»,ROOTR<4),CURNT(3» 
DIMENSION RM3t2)»ALPHAA(2),THETAA(3,2).0MEGAA(2) 
DIMENSION R3(3*2)*ALPHAB(2)*THETAB(3.2) tOMEGAB(2) 
COMPLEX R00TAI4) *CKA(3*4ltCNUMA(3,4) ,DN0MA<4I ,S 
COMPLEX R00TB(4>.CKe<3f4l*CNUMB<3,4l#DN0MB<4> 

COMP LEX CMPL X 

REAL LltNTI*L2INTI»Ll,L2,LlS.KUK2,LCl*IDUMP 
INTEGER CYCLE 


COMPONENT VALUES. 

Cl * 2 640. OE- 06 
C2 > 880.0E-06 
C3 * 180. E-6 
C6 « 1.1E-9 
LI » 48.0E-06 
L2 * 480.0E-06 
LC1 * 0.003E-06 
R1 * 85. OE 03 
R2 ■ 4.99E 03 
R3 - 15.22E 03 
R12 - 750.0E 03 
R20 » 50. OE 03 
R21 » 400. OE 03 
R22 * 7.84E 03 
RD = 1.0E-03 
RL - 4.00 
RL1 - 10.0E-03 
RL2 * 60.0E-03 
RC1 - 10.0E-03 
RC2 * 10.0E-03 

INITIAL VOLTAGES AND CURRENTS. 

VD = 0.6 
VIN * 60.00 
ClINTV » 29.50 
C2INTV * 29.50 
L1INTI * 0.0 

C THE FOLLOWING IS AN APPROXIMATION FOR L2INTI AND C3INTV. 
L2INTI * C2INTV/RL 
C3INTV - 0.90*R2/ 1 R1+R2 1 *L 2 1 NT I 
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C 


c 


c 

c 

c 


c 


c 


VC1 « C1INTV 
VCURNT * C3INTV 
CURNT2 a ABSU2INTII 

CYCLE » 0 
MAXC - 500 
N - 1 
T * 0.0 
TTOTAL * 0.0 
K1 » 0.70E-06 
K2 * 1.08 


WRIT £(6. 3801 VIN 
WR IT £( 6.390 1 PL 

DETERMINE ROOTS OF FILTER* 


XCOF ( 1 » 
XCOF ( 2 ) 

XCOF ( 3 I 


RL+RL2+RL1+RD 

Ll+L2+Cl*( RCl*RL2+l RLl+RDl*RCl+( RL 1+RD) *RL2 
+RL*RCl*RL*(RLl+ROmC2*(RL*RC2+RL*RL2+RC2* 
RL2*RL*(RL1«-R0I+RC2*(RLU-RC» ) 

C2*L2*( RL+RC2 I*C2*L1*( RL*RC2 )+C 1*L1* ( RL+RC1 
+RL2»*Cl*L2*(RCl+RLl+R0>+Ci*C2*<RL*RCl* 
RC2*RL*RL2*RC1*RC1*RC2*RL2*-RL*RC1*(RL1+RD> 
♦RL*R;2*|RLl4RDH-RL*RL2*(RLl+RDI + RCt*RC2* 


IRLH-RD)+RC2*RL2*(RL1+R0) 1 
XCOF (41 » C1*L1*L2*C1*C2*L1*(RL*RC1+RL*RC2+RL*RL2+ 

RC1*RC2«*RC2*RL2M-C1*C2*L2*(RL*RCI+RL*(RL1+ 
R0I+RC1*RC2+RC2*(RL1+RD) ) 

XCOF (5 1 » Cl*C2*Ll*L2*(RL«-RC2> 

XCOF A l « XCOF ( 1 I 
XC0FA2 * XCOF ( 2 ) 

XCOF A3 » XCOF ( 3 ) 

XC0FA4 * XCOF ( 4 ) 

XC0FA5 * XC0F<5> 


WR IT E ( 6f 400) 

WRITE ( 6*410 ) 

WRITE(6,420> 

H * 4 

CALL POLRT ( XCOF . COF* M, ROOTR * ROOTI * I ER I 
DO 100 1-1,4 

ROOT A(I) * CMPLXt ROOTR ( I ) , ROOT I (HI 
WRIT E( 6*4301 I , ROOTR (I> .ROOTI (I> 

100 CONTINUE 

00 130 I* 1 i 4 
S * ROOT A ( I ) 
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DN04AI I) * XCOFI 51 
K = 0 

liO K = K + l 

I F ( K *EQ. 1 I GO TO 110 
IF< K .GT. A- I GO TO 120 
ONOMA<1> - DNQMA ( l ) * ( S-ROOT AIK)) 

GO TO 110 
120 CONTINUE 
130 CONTINUE 

XCOFUI * 1*0 

XCOFI21 * C2*RL+C2*RC2*C1*RL+CI*RC1+C1*RL2 
XCOF (3) * C1*L2+C1*C2*RC1*RL+C1*C2*RC1*RC2*C1+C2*RL2* 
6 RL+C1*C2*RL2*RC2+C1*C2*RC2*RL 

XCOF (4) * C1*C2*L2*RL+C1*C2*L2*RC2 
XCOF I 5 1 * 0*0 
C 

WR IT E ( 6. 400 1 
WRITE ( 6.440 ) 

WRITEI6.420) 

C 

H a 3 

CALL POLRTIXCOF.COF.M.ROOTR.ROOTI ,IER) 

00 140 1*1*3 

ROOT B( I ) = CHPLXI ROOTRI II* ROOT KIM 
WR I T E < 6* 430 1 I * ROOTR (1)* ROOTI I I i 
140 CONTINUE 

00 170 1*1.3 
S - ROOTBII) 

ONOMBII) * XC0FI41 

K - 0 

150 K * K+l 

IF( K .EQ. I ) GO TO 150 
IF! K *GT. 3 I GO TO 160 
ONOMBII) * DNOHBI I )* I S-ROOTBt K ) ) 

GO TO 150 
160 CONTINUE 
170 CONTINUE 
C 

WR IT E (6.450) 

WRITE I 6.4601 C 1INTV. L II NTI 
WRITEI6.470) C2INTV.L2I NTI 
WRIT E( 6.480 ) C3INTV 
C 

160 CONTINUE 

C WITH STATEMENT 180. THE REGULATOR CYCLE BEGINS. 

C 

C CYCLE EQUALS THE CYCLE NUMBER THE REGULATOR IS ABOUT 
C TO GO THROUGH. 
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C 

CYCLE * CYCLE+1 
WRITE! 6*490 I CYCLE»TTOTAL 
C 

00 190 1*1,3 
S * R00TBI I) 

190 CNUMBI2,!) * S**2*I C1#C2* I RL+RC2 ) *L2*L2 INTI )+S* I C 1* 

C L2*L2!NTI+C1*C2*( I RL+RC2 ) *C1 INTV-Rl* 

& C 2 INTV I ) ♦C 1*C l INTV 

00 200 I * 1,3 
S - R00TBI I) 

200 CNUi 8! 3, 1 1 - S**2*!C1*C2*IL2*C2INTV+RC2*L2*L2INTI ) ) + S# 
t ! C 1*C2*( < RC 1*RL2I *C21NTV+RC2*C 1 INTV ) +C 1 

£ *L2*L2INTI >MCI*C1INTV+C2*C2INTV) 

NOW CALCULATE THE COEFICIENTS, FREQUENCIES ANO 
DAMPING FACTORS OF THE LDOP CURRENTS WHILE THE 
FILTER IS IN MODE B. 

00 240 I * 2,3 
00 230 J * 1,3 

IF( A IMAGI ROOTB! J ) I ) 230,220,210 
210 CKB(I,1) * CNUMBI I, J t/ONQMB! J) 

RBI I , 1 ) - CABSI CKBI 1,11) 

TOPS * AIMAGI CKBI 1,11) 

BOTB * REAL I CKBI I , 1 ) ) 

THETABI 1,1) * ATAN2I TOPB,BOTB) 

OWES ABU) * AIMAGI ROOTB I J I ) 

ALPHABU) * -REAUROOTBOn 
GO TO 230 

220 CKBI 1,2) * CNUMBt I , J ) /ONOMBI J ) 

RBI 1,2) * CKBI 1,2) 

ALPHABU) * -ROOTB! J ) 

230 CONTINUE 
240 CONTINUE 

NOW CALCULATE T DWELL* 

TOWELL = 1 1 VCURNTM T*5E3*R3)-15*0*R3) /7*5E3*VC1 ) 

L +C1/CURNT2 

IF! TDWELL .LT. 0.0 ) TDWELL-0.0 
T * TOWELL 

NOW CALCULATE THE FILTER CURRENTS ANO VOLTAGES 
DURING TOWELL. 

HH * THETABI 2, l l-AT AN2I0MEGAB I 1 > * -ALPHAS I 1 ) I 
SQB = SQRT I ALPHABI 1 ) **2+0MEGAB I 11**2) 
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C 


C 


OTl * T/5.0 
CURVTI 1) * 0.0 
00 250 K* 1,5 
T » FLOAT ( K | *0T1 


CHARGI * 0.0 

CUPVTI 2) = 2.*RB(2«L l*EXP(-ALPHAB ( 1 ) *T ) *C0S ( OMEGAS { 1 ) 
£ *T*THETAB<2#in*RB(2,2)*EXP(-ALPHAB(2>*T> 

CURNT ( 3 ) * 2**RB( 3, 1 1*EXP(-ALPHAB< 1 l*T )*CQ$< OMEGABI I) 
£ *T*THETAB<3,1) >*R8<3,2)*EXP(-ALPHAB(2)*T) 

VC1 * -l/Cl*(2*RB<2,l)*EXP<-ALPHABm*T|/$QB* 

& COStOMEGABt U*T+HH> +RB (2,2 Ui -ALPHAS < 2 M * 

£ £XP(-ALPHABI2I*T| t+l.OSC l* < 2.0*RB ( 2 , 1 ) / 

£ $QB*C0S(HH)*RB<2,2>/ (-ALPHAS! 21 ) >*C1INTV- 

£ RCL*CURNT<2) 

VOUT = RL*CURNT ( 3 ) 

TTOTAL * TT0TAL+0T1 
N * N + l 
250 CONTINUE 


WRIT E 1 6* 5001 
WRIT E ( 6*510) 
WR ITE ( 6. 520 ) 
WR IT E (6, 530 ) 


TOWEL L 

CURNT(l) .CHARGI 
CURNTf 2) * VC1 
CURNTI 3) , VOUT 


NOW CALCULATE TCHARG. 

VI * 0.20*CURNT2 

THIVG » ((R12+R20)/(R2l+R22)*(Vl+0. 6)+VI N+0. 6 ) / 1 VI N-Vl ) 
TCHARG « ( R12+R20 )*( R21+-R22) *C6/( R12+R20*R2l*R22 ) 

£ *ALOG( THING) 

NOW CALCULATE THE FILTER CURRENTS AND VOLTAGES 
DURING TCHARG. 

CHARGI « 0.0 
T = TCHARG 
TOLD * 0.0 
T 1 * 0.0 
K « 1 

0T2 * T/50.0 
00 260 I s l» 50 

US » 4.0*<L1-KI*CHARGI**K2) 

DELI ■ V IN*DT 2/LIS 
CHARGI * CHARGl+DELl 
T 1 * T1+DT2 

I F( I .EQ. 1 ) GO TO 255 
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IF! I *NE. 10*K | GO TO 260 
K * K* 1 
255 CONTINUE 

T * TDWELL*T I 

CURNT!2» » 2«*RB( 2« l )*EXPt -ALPHAS! 1 1 *T I *COS( OMEGA B ( 1 1 
£ *T+THETAB( 2* I ) »*RB!2.2I*EXP(-ALPHAB!2»**» 

CURVT(3> * 2.*RB!3.l»*EXP(-ALPHAB! 1»*T»*COS!OMEGAB!1) 
6 *T+THETA8!3.1I **R8< 3. 2 I *EXP!-ALPHAB( 2 1 *T ) 

VCI * -1/C1*(2*RB(2.1»*EXP(-ALPHAB(1>*TI/SQB* 

£ COS (ONE GAB! ll*T*HH (+RB! 2 .21 /(-ALPHA8(2) I * 

£ EXP (-ALPHAS! 2 »*T> ) *1 *0/C 1* ! 2* 0*RR! 2 , U / 

£ SQB*COS(HH) *RB! 2.2) /!-ALPHAB!2) ) ) *C1INTV- 

£ RCl*CURNT(2t 

VOUT * RL*CURNT( 31 
TTOTAL » TT0TAL+T1-T0LD 
TOLD « Tl 
N * N*1 
260 CONTINUE 

CURST2 = ABSI CURNT 12)) 

VCURNT * 0*90*CURNT2*R2/(R1+R2»-!0*90*CURNT2*R2/(R1 
£ *R2)-C3INTV>*EXP(-(Rl*R2>/(Rl*R2*C3l*T> 

C 

WRITE!6.540) TCHARG 
WRITE(6,510> CURNT ( 1 ) .CHARGI 
WR IT E ! 6. 520 I CURNT ! 2 ) .VC1 
WRITE! 6*530 ) CURNT! 3 1 . VOUT 
C 

C UPDATE INITIAL CONDITIONS* 

C 

L1IMTI * 4*0*CHARGI 
L2IVTI * CURNT! 21 
C1IVTV * VCI+RC1*CURNT(2I 
C2INTV « VOUT 
C3IVTV * VCURNT 
CURNT1 a CHARG! 

CHARGI » 0*0 
C 

C NOW THE INPUT INDUCTOR HAS BEEN CHARGED AND IS GOING TO 
C DUMP ANO THE FILTER WILL BE IN MODE A. 

C 

DO 270 1-1.4 
S « ROOT A! I l 

270 CNUMAII.I) » s**3*( C 1 *C 2 *L 2 *( RL*RC2)*L 1 *L l INTI I +S** 2 * ! 
£ C1*L2*L1*L1INTI*C1*C2*(RL*RC1*RL*RC2*RL* 

£ RL2+RC 1*RC2+RC2*RL2I*L1*L1 INTI -C1*C2*L2* 

£ !RL+RC2>*CIINTV-C1*C2*L2*(RL*RC2»*VD*C1* 

£ C2*( RL*RC1*RC1*RC2I*L2*L2INTI »*S*(C2*!RL 

£ +RC2I*H*HINTI+C1*RL*L1*L1INTI*C1*(RC1* 
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C 

C 

C 

C 


£ 

£ 

£ 

£ 

£ 

£ 

& 

£ 

& 


Rl2J*LI*LnNTI-Ci*L2*eiINTV-Cl*C2*(RL*ftC2 
+RL*RL2*RC2*RL2»*CXINTV-CX*L2*V0-C1*C2*< 
Rt*RC l+RL*RC2+RL*RL2+RC X*RC2*RC 2*RL2) *VD 
+C2*(RL+RC2I*L2*L2INTI+CL*RCI*L2*L2INTI- 
C1*C2*RC1*RL*C2INTV) ♦ (L1*LXINTI+L2* 
L2INTI-C X*(RL*RL2i*CX INTV*C2* (RL*RC2) *VD 
-CX*(RL+RCU-RL2)*VD-C2*RL*C2INTV) 
-<-V0/IRL+RL2+RLX*RD) )*($**3* 
XCOFA5+$**2*XCOFA4+S*XCOFA3+XCOFA2) 


00 280 1=1,4 
S * ROOTAt I > 

280 CNUWAI 2, 1 1 * S**3*tCX*C2*LX*tRL+RC2»*L2*L2INTI >+$**2* 

& <CL*Li*L2*L2INTM-CX*C2*LX*<RL+RC2>*ClINTV 

& -C1*C2*L X*RL*C2INTV+C X*C2* (RL*RCX*RL* (RLX 

6 ♦RD>+RCX*RC2+RC2MRLX+RD) »*L2*L2INTI+C1* 

£ C 2 *RCl*<RL*RC 2 )m*LlINTn+S*IC 2 MRL+RC 2 

6 »*L2*L2INTI+CX*tRCI*RLI+RD>*L2*L2!NT!+CX 

£ *C2*(RL* (RL1*RD)+RC2*( RLX+RD) >*CIINTV-CL 

6 *C2*(RL*RC1+RL*(RL1+RD* )*C2INTV+C2*(RL*RC2 

£ >*LmiINTI+Cl*RCim*LlINTI-Cl*C2*RCl*< 

£ RL*RC2 l*VD+Cl*Ll*C!INTV>+( LX*LLINTI+L2* 

£ L2INTI-C2*RL*C2!NTV+C l*IRLi+RDI*Cl I NTV-C2 

£ *<RL«-RC2l*VD-CX*RCl*VD» 

£ -<-VD/«RLX+R0+RL2*RU >*<S**3* 

£ XC0FA5*S*#2*XC0FA4+S*XC0FA3+XCQF A2 > 

00 290 1-1*4 
S * ROOT A ( 1 1 

290 CNUNA(3*I» * S**3*CCI*C2*LX*L2*C2INTV+CI*C2*LI*RC2*L2 


£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 


♦L2INTI I ♦S**2*(CX*C2*HM RCX+RL2 I *C2I NTV 
♦Cl*LI*L2*t2INTt+CX*C2*H*RC2*CiINTV+CX*C2* 
L2*« RCX+RLL+RD)*C2INTV+C1*C2*<RCI+RLM-RD» 
*RC2*L2*L2INTI*C1*C2*RCX*RC2*LX*LI INTI > + 

$*( C2*L2*C2INTV+C2*RC2*L2*L21NTI+C2*L X* 
C2INTV+CXnX*CXINTV+CI*C2*<RCI*RL2+RCX*<RLl 
+R0) +RL2*< RLX+RD) ) *C2INTV+CX* ( RCX+RLl+RD) * 
L2*L2INTH-CX*C2*RC2*(RLI+RD>*CXINTV+C2*RC2* 
L1*H INTI +CX*RCX*LX*LX INTI -CX*C2*RCX*PC2* 
VD) + (LX*HINTI+L2*L2INTI+C2*RL2*C2INTV+C2 
*<RH+R0)*C2INTV+CX*<RLX + RD>*CXINTV-C2* 
RC2*VD-CX*RCX*VD) 

-(-VD/(RLX*RD+RL2+RLI )*<S**3*XC0FA5+S**2* 
XC0FA4+S*XC0FA3+XC0FA2» 


NOW CALCULATE COEFICIENTS# FREQUENCIES AND DAMPING FACTOR 
OF THE LOOP CURRENTS WHILE THE FILTER IS IN MODE A. 


DO 340 1*1*3 
K * X 
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00 330 J*l,4 

IF< AIMAG(ROOTAIJI) ) 330*330,300 
300 GO TO (310*320) ,K 
310 CKA< t * 1 1 * CNUMA( I*JI/DN0MA( J) 

RA(I * 1 1 * CABSICKAt 1*11) 

TOPA * AIMAGt CKA( I, II I 
BOTA « REAL (CKA( 1*1) ) 

THET AA( 1 , 1 ) * ATAN2( TOPA* BOTA) 

K * 2 

IFU .GT. It GO TO 330 
OMEGAA(l) * A I HAG (ROOTA ( J ) ) 

ALP-tAA(l) = -REALIROOTA* J) ) 

GO TO 330 

320 CKA( 1*2) * CNUMAC I,J)/DNOMA<JI 
RA(I ,2) * CABS(CKA( 1*2) t 
TOPA * AIMAGCCKAI 1,2) ) 

BOTA * REALICKA ( 1*2) ) 

THET AA( 1,2) - ATAN2IT0PA,B0TA) 

IFU .GT • II GO TO 330 
OMEGAAI 2) = AIMAGIROOTAI J ) ) 

AL PHAAI2) * -REAL(ROOTAU) I 
330 CONTINUE 
340 CONTINUE 

NOW CALCULATE TDUMP. 

K * l 

T2 * 0.0 

DT3 = 0.0 

TOLD * 0.0 

DELI * CURNT1/40.0 

DO 360 1*1,41 

LIS * 4.0*(U-K1*<DEU*(41-IH**K2) 

T2 * T2+DT3 

DT3 * 0EL1*L1S/I VC1+0.6I 
I F ( 1-1 .LE. 2 I GO TO 350 
IF( (1-1) .NE. K*5 I GO TO 360 
K * K*l 
350 T * T2 
TDUMP * T 

NOW CALCULATE FILTER CURRENTS AND VOLTAGES MHEN THE INPUT 
INDUCTOR BEGINS TO DUMP. 

CURNTUI * 4.0*1 DEL I * 1 4 1 — I 1 1 
IDUMP » CURNT ( 1 1 

CURNTUI * 2.*RA( 2* 1 )*EXP(-ALPHAA( 1 )*T l*COS( OMEGAA ( 1 ) *T+ 
£ THETAAI2, 1) l+2.*RA(2,2)*£XPI -ALPHAAI 2 )*T)*C05 



o o o n o o o o r> 


65 


£ t OMEGA A < 2I*T<-THETAAt2,2) > -VD/t RL1 *RQ *RL2*RL I 

CURMT(3I = 2.*RM3»n*EXP(-ALPHAA<n*T)*COSO*EGAAtn*T* 
6 THETAAi3.ll ) +2. *RA( 3 ,2 > *EXP t -ALPHA A t 2 1 *T > *COS 

6 tOMEGAAi2)*T+THETAA(3,2n-VD/<RLl + RO + RL2+RU 

VI Cl * -LC1*4.0*DELI/0T3 

IF( T *LT • 2.0E-06 ) VLC l»LC l*Ll I NT 1 / 1 . 0E-06 
VC1 = 1.0/Cl*<-4.0*DELI/t2.0*DT3i*T**2MLUNTI- 
£ CURNT (21) *Tl +C 1 INTV+RC l* (CURNTi 1 > -CURNT 1 2 ) ) +VLC 1 

VOUT « RL*CURNT <3 ) 

TTOTAL - TTOTAL+T-TOLD 
TOLO = T 
N > N+l 
360 CONTINUE 
C 

CURNT2 - ABSiCURNT (2)1 

VCURNT » 0.90*CURNT2*R2/(R1+R2)-<0.90*CURNT2*R2/(R1 
£ +R2>-C3INTV)*EXPi-IRl+R2>/<Rl*R2*C3>*T) 

HR IT E ( 6.550 ) TDUMP 
WRI T E < 6, 510 ) CURNTi 1 ) .CHARGI 
WRIT E (6.520 1 CURNT(2),VC1 
WR IT E 16.530 I CURNTi 3 )» VOUT 

UPDATE INITIAL CONDITIONS. 

LI INTI * CURNTi 1 ) 

L21MTI * CURNTi 21 
Cl INTV * VC1+RCI*CURNT(2> 

C2INTV * VOUT 
C3INTV * VC UR NT 


THE CHARGING INDUCTOR HAS DUMPED AND THE LATEST FILTER 
CURRENTS AND VOLTAGES HAVE BEEN CALCULATED. NOH RETURN 
TO STATEMENT 180 AND START ANOTHER CYCLE OR END. 

I F ( CYCLE .GE. MAXC ) GO TO 370 
GO TO 1B0 
370 CONTINUE 

380 FORM AT ( *2*.25X,* INPUT VOLTAGE - VIN * *,F1U4) 

390 FORM AT i *0* . 25X, 'LOAD RESISTANCE = RL * *,F10.4> 

400 FORMAT! *-•*///» 17X. 1 THE FOLLOWING ARE THE ROOTS TO* . 

£ • THE OUTPUT FILTER*) 

410 FORMATi* *.16X.*WHILE IT IS IN MODE A.*) 

420 FORMATi *-• ♦ 34X. * REAL PART * . 7X, * IMAG PART* ) 

430 FORMATi *0' ,23X, *ROOT(*. II, *) * * , E12, 5 ,4X , E12. 5 ) 

440 FORMATi* '.UX^WHILE IT IS IN MOOE B. ' ) 
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450 FORMAT ( ' 1*»18X* 'STARTING INITIAL CONDITIONS'! 

460 FORMAT! •-* ,21X* 'C1INTV * * , 6 l 1.4* 4X, • L II NTI = ',£11.41 
470 FORMAT! * ' , 21 X, ' C2 INTV « ■ *E1 1. 4, 4X, * 1 21 NT I = »,FU.4) 
480 FORMAT ( * ' »21X* 'C3INTV * ' # EU.4,///> 

490 FORMAT ('0'»15X»' BEGINNING CYCLE NUMBER ',I4,5X, 

& 'TTOTAL « SE11.4) 

500 FORMAT! '0» t 21X, 'VALUES AT THE END OF TDWELL « '.Ell.41 
510 FORMAT I * • ,2'SX* 1 CURNT (X ) * • .EU.4.2X, 'CHARSI = »»E11.4i 
520 FORMAT! * • ,25X» *CURNT(2) * * t Ell. 4, 5X f • VC1 = •,EU.4» 

530 FORMAT ! • • ,25X, 'CURNTI3I * * .E11.4.4X, *VOUT « ••EU.4) 
540 F0RMAT('0*,21X, 'VALUES AT THE ENO OF TCHARG » «,E11.4» 
550 FORMAT !*0*»2XX*' VALUES AT THE ENO OF TOUMP = ',E11.4) 

STOP 

ENO 
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INPUT VOLTAGE * VIN = 
LOAD RESISTANCE = Rl = 


THE FOLLOWING ARE THE ROOTS TO THE 
WHILE IT IS IN MODE A. 

REAL PART 

ROOTtll « -0.21639E 03 
ROOT (21 * -0.21639E 03 
R00T(3) * -0.22736E 03 
R00TI4) # -0.227366 03 

THE FOLLOWING ARE THE ROOTS TO THE 
WHILE IT IS IN MODE 8. 

REAL PART 

ROOT ( 1 ) * -0.706536 02 
ROOT ( 2 1 * -0.18967E 03 


60.0000 

4.0000 

OUTPUT FILTER 

I MAG PART 
0.143856 04 
0.L4385E 04 
0.29852E 04 
0.29852E 04 

OUTPUT FILTER 

I MAG PART 
0.0 

0.17690E 04 


R00TI3I * -0. 1 89676 03 


0.17690E 04 
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STARTING INITIAL CONDITIONS 


C 1 1 MTV = 0.2950E 02 L1INTI * 0.0 

C2INTV * 0.2950E 02 L2INTI * 0.7375E 01 

C3INTV = 0.368 IE 00 


BEGINNING CYCLE NUMBER 1 TTOTAL * 0.0 

VALUES AT THE END OF TDWELL * 0.6263E-04 

CURNTU) * 0.0 CHARGI * 0.0 

CURNT ( 2 ) = 0.72976 01 VCX » 0.2925E 02 

CURNT (3) » 0.73746 01 VOUT * 0.29506 02 

VALUES AT THE END OF TCHARG = 0.2967E-04 

CURNTU) * 0.0 CHARGI * 0.1008E 02 

CURNT ( 2 ) * 0.7252E 01 VC1 * 0.2917E 02 

CURNT ( 3) = 0.7373E 01 VOUT » 0.29496 02 

VALUES AT THE ENO OF TDUMP = 0.5872E-04 

CURNTtl) * 0.0 CHARGI = 0.0 

CURNTU ) * 0.72I5E 01 VC 1 * 0.2950E 02 

CURNT (3 ) * 0.7371E 01 VOUT * 0.2948E 02 

BEGINNING CYCLE NUMBER 2 TTOTAL = 0.15106-03 

VALUES AT THE END OF TDWELL » 0.6337E-04 

CURNTtl).® 0.0 CHARGI = 0.0 

CURNTU ) * 0.7150E 01 VC1 * 0.2933E 0? 

CURNTU) * 0.7367E 01 VOUT = 0.2947E 02 

VALUES AT THE ENO OF TCHARG * 0.2922E-04 

CURNTtl) = 0.0 CHARGI * 0.9912E 01 

CURNT { 2) = 0.7U3E 01 VC1 = 0.2925E 02 

CURNTU) * 0.7365E 01 VOUT = 0.2946E 02 

VALUES AT THE ENO OF TOUMP * 0. 5771E-04 

CURNTU) * 0.0 CHARGI = 0.0 

CURNTU) = 0.7090E 01 VC 1 * 0.2956E 02 

CURNT 1 3 ) * 0.73606 01 VOUT * 0.2944E 02 

BEGINNING CYCLE NUMBER 3 TTOTAL ® 0.3013E-03 



V. THE AC OUTPUT IMPEDANCE 


It is possible to calculate the ac output impedance of a regu- 
lator as a function of frequency and use the result to obtain an insight 
into the regulation of the system. The regulator's output stage is a 
low-pass filter usually composed of one or more rc-filter sections each 
of which is composed of inductors and capacitors having soaie series 
equivalent resistance. The output impedance will exhibit resonances , 
or zeros and poles of impedances , at certain frequencies* At the poles , 
the impedance of the filter becomes relatively large and may reach a 
value that will lead to poor voltage regulation. Thus , a design require- 
ment is to choose component values that will place the resonance points 
outside the range of any frequencies that might be introduced by fluctu- 
ating loads . 

One common method of calculating the ac output impedance is to 
inactivate all independent sources and mathematically apply a test 
voltage, v T (f,t), at the output. The scripts (f,t) indicate that the 
voltage v,j,(f,t) will be a function of frequency and time. The ac out- 
put impedance is equal to the ratio of the phasor equivalents of the 
voltage applied at the output and the resulting current. When linear 
feedback signals are present in the regulator, they must be accounted 
for by realizing that they are dependent sources and must be included 
in calculating the output impedance which can be done by using linear 
analysis techniques. 
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When nonlinearities are present, analysis becomes more compli- 
cated and often leads to tedious calculations which force one to seek 
appropriate simplifications . Consider now the CBRM regulator which has 
two sources of nonlinearity as shown in Figure 5-1, a simplified model 
of the regulator. The nonlinearities are due to the switching modulator 
and the circuitry that generates a feedback voltage, FB 2 , proportional 
to ' the absolute magnitude of the current in the second loop of the out- 
put filter. Xn this system then, an unrestricted linear analysis can- 
not be used due to the presence of frequencies, other than the original, 
produced by these nonlinear devices. Rather than perform a nonlinear 
analysis on the CBRM regulator, a piecewise-linear analysis will be done. 
To justify this, it will be shown that the influence of the nonlinear 
feedback signal, FB 2 , becomes negligible when the frequency of v^,(f,t) 
exceeds a low-frequency range, and it will also be shown that the switch- 
ing modulator can be linearly modeled. In this manner, a restricted 
model will be produced for which the ac output impedance will be develop- 
ed. Before covering the nonlinear feedback path, the switching modu- 
lator will be discussed. 

Since the switching frequency is much higher than the cutoff 
frequency of the output filter, it is assumed that the ripple components 
reaching the feedback connections are small and any ripple components 
reaching the input to the modulator can be neglected. With this assump- 
tion, the output voltage of the switching modulator can be replaced by 
its average value as far as its effect on the output and feedback volt- 
age is concerned. As discussed in Chapter II, the average value of the 
switching modulator output voltage is a function of its duty cycle which 




Figure 5-1. A simplified nonlinear model of the regulator. 



72 


in turn is a function of the modulator's input voltage. Based on the 
above discussion, the switching modulator can be modeled as a voltage- 
controlled voltage source to obtain the linear model shown in Figure 
5-2, The Justification for eliminating the nonlinear feedback path will 
now be given. 

- Considering the nonlinear feedback path in Figure 5-1, the 
feedback signal, FB 2 , is a voltage resulting from monitoring the magni- 
tude of the current in the second loop of the filter. The path gain, A 3 , 
has a value of approximately 0.55 and a cutoff frequency of approximately 
0.2 Hz. The feedback signal FB^ is a voltage proportional to the voltage 
v c and is summed with a fraction of the reference voltage. The gains 
A^ and A^ are produced by a voltage divider network and have values of 
approximately 0.67 and 0.33, respectively. 

For frequencies greater than ten times the cutoff frequency , 

0.2 Hz, the feedback signal FB^ becomes small and its influence on the 
system is negligible when compared with FB^, Since v c will be large at 
low frequencies, FB^ will be large. At higher frequencies, FB^ will 
decrease but so will FB^ • Thus, FB^ will have the dominant influence 
once the cutoff frequency of the FB 2 path is exceeded. However, note 
that there is a high frequency range reached for which all feedback will 
be negligible and the output impedance will be determined by the resist- 
ance and inductance of the capacitor in the filter's output branch. 

This is true since the capacitors offer little Impedance and the induct- 
ance of L2, a filter inductor, is many times greater than which is 
the small series equivalent inductance of C2. If the frequency of 

is restricted to values greater than the 0.2 Hz cutoff value, 
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the nonlinear feedback path can be neglected leaving the linear system 
shewn in Figure 5-2. 

Before the ac output impedance can be calculated for the restrict- 
ed linear model * the value of must be determined. To do this* it is 
recalled from Chapter IX that the purpose of the current monitor and 
associated circuitry was to force a reduction in output voltage as a 
function of output current. The resulting characteristic is shown in 
Figure 5-3. This response* good for slow changes in load current* yields 
an effective dc resistance of 213 milliohms. Using this value and solv- 
ing ratio Vip(0,t)/i>].(0*t) produces a value of approximately 100 for A^. 

Now solving for the ac output impedance* Z(w) , of the circuit 
in Figure 5-2 yields: 


ZCw) 


B lV B 2 + B 3 ) + B 2 B 3 B 4 <1 + K) 

+ B 3 + B*) + B 2 (B 3 +”B 4 )"C1 + K) 


(5-1) 


where K = 

= .Rl-l + jwLl 

B 2 " R C1 + j(tl)l, Cl " 1/toC1) 

B 3 “ ^2 + jwL2 

B 4 « ^ + j (a)L c2 - l/wC2) , 


and 
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OUTPUT CURRENT (AMPS) 


Figure 5-3. The output voltage/output current characteristic showing 
the programmed slope. 
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The magnitude of the ac impedance was calculated and plotted 
for a range of frequencies and is shown in Figure 5*4. Also shown on 
the same plot is the magnitude of the impedance due just to the output 
filter* i*e. * with no feedback. 







VI. POWER SHARING 


Now that the single regulator has been modeled and analyzed, a 
study will be made of the complete power system in which eighteen such 
regulators are paralleled and are under the control of a power sharing 
feedback network. 

It was found in the previous analysis that the Thevenin voltage 
of the regulator was approximately 30. A volts and the Thevenin resist- 
ance was approximately 213 milllohms. Thus, the general expression for 
the output voltage, in this chapter called V^, as a function of the load 
current, 1^, is: 

V 1 * V K2 - Vm < 6 -D 

where V ^ and R & is the Thevenin voltage and Thevenin resistance, respect- 
ively* The value of is in the range of 27 to 30 volts, depending on 
the current demand which could be as large as 15.5 amperes. 

Figure 6-1 illustrates how one regulator is connected to the out~ 
put power bus and is under the influence of the power sharing circuit* 

With no feedback, i.e., Vj equal to zero, the regulator portion of the 
figure is a model for Equation 6-1. The feedback voltage, V^, is devel- 
oped as follows. The output bus voltage, Vg, is compared with a refer- 
ence voltage, V R1> of 29.5 volts and the difference is amplified by 
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which has a value of 10.8 amperes per volt. The output of is regard- 
ed as a programmed current, Ip, and is used as a basis for comparison 
with the actual measured output current, I m . The difference in I and 

r 

Ijjj is amplified by A£» whose value is 3.35 volts per ampere, and the 
resultant is Vf. Thus : 

V f - A^A^Vm - V B ) - I m ] . (6-2) 

The regulator output voltage then, with the feedback voltage, Vf, is 

V 1 “ V R2 " Vm + V f • < 6 " 3 > 

The output bus voltage, V will be equal to minus the volt- 
age drop across R^, the distribution line resistance. 

Any number of regulators can be paralleled as shown in Figure 

6-2. Considering the network, Vf is the output of the regulator, 

tin 

R U is the distribution line resistance of this i network, and 1^ is 
the current this network furnishes to the load. Each network has as amp- 
lifier with a gain of which will be considered the same for all regu- 
lators. To insure the same programmed current for all regulators, one 
reference voltage, V Rl ' and one amplifier, Af, are used. The sum of all 
the I m currents equals Xy, the total current delivered to R^, the load 
resistance. In the Skylab/ATM power system, there are eighteen regulators 
in parallel, but to conserve the generality of the following analysis, it 
will be assumed that n such regulators are in parallel. 
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A more useful feedback diagram will be constructed based on a 
relation for the bus voltage of n regulators in parallel, and the rel- 
ationship obtained from this diagram will reveal the interaction among 
the regulators. The total output current, I^, can be written as: 


■v, 


M 


ji _ I 


®L 


ml 


+ I 


m2 


+ I 


mi 


+ I 


an 


( 6 - 4 ) 


The individual currents, such as , can be replaced by an equivalent 
form. 


I 


mi 



( 6 - 5 ) 


Thus: 



(6-6) 


Solving for Vg yields: 




Z v i/%i 

1 + ^ Ji 1/Rii 


( 6 - 7 ) 
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Referring to Figure 6-2, it is seen that: 

V i " v R2i + A 1 A 2 v R1 ~ a 1 a 2 V B ' (Bsi + a 2> ^mL • (6-8) 

Using Equation 6-5 and 6-8 in 6-7 yields: 


V B 


p £ v R2i + a 1 a 2 V R1 

^ Ji »U * Si + *2 

1 + S. £ v 1 + . A i A : . 

^ i»l %i + *si + A 2 


(6-9) 


By substituting Equation 6-8 into 6-5 » an expression for the 
individual regulator output currents is obtained. 


v R2i + A 1 A 2 V R1 " (1 + a 1 a 2) v B 
hi + hi + a 2 


( 6 - 10 ) 


Using Equation 6-9 , the final feedback diagram for n regulators 
in parallel is constructed and shown in Figure 6-3. The individual regu- 
lator output currents as well as the total output current are labeled. 

To observe the relationship between bus voltage and output cur- 
rent, Equation 6-9 was used to calculate for the possible range of 

load resistance. Identical regulators were assumed with line resist- 
ances of 50 mil li ohms. The results for different numbers of regulators 
in parallel are plotted in Figure 6-4. Noting the slope of the V B /I M 




Figure 6-3. The final feedback diagram 



r n regulators in parallel. 
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Figure 6-4. Regulation of paralleled regulators under 
power share control. 


300.00 



86 


curve for eighteen in parallel, the system Thevenin resistance is approx- 
imately 5.6 milliohms. Thus, the new Thevenin resistance for each reg- 
ulator, to be called is 18 times 5.6 milliohms or 100 milliohms. 

s 

This can be shown mathematically by rewriting Equation 6-10 as: 

v 1 A 1 A 2 „ R U + ^i + A 2 . ,, . 

V B = 1 + AjA 2 V E2i + 1 + AjA 2 Vr 1 " 1 + A^A 2 * (6 “ 11) 


How taking the derivative of Vg with respect to I mi . yields the new 
Thevenin resistance, R^, of: 


R 


‘si 


%i + ^i + a 2 
1 + A^ 


* 97 milliohms. 


(6-12) 


Figure 6-5 shows how the current in one regulator of 18 in par- 
allel varies as a function of its line resistance. To create this fig- 
ure, the line resistances of 17 of the regulators were set to 50 tnilli- 
ohms and the resistance of the eighteenth was varied from 10 milliohms 


to 110 milliohms. 
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(b) r £18 (milliobms) 


Figure 6-5. (a) The output current of each regulator and the 

total output current, Iy. (b) An expanded view 
of the individual regulator currents. 
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VII. RESULTS AND CONCLUSIONS 


This chapter presents some of the results that were obtalmed 
through the computer program model of the regulator. The results are in 
the form of graphs created from data generated by the program and are an 
indication of the type of information that can be obtained. 

First, for a variety of load resistances, the different switch- 
ing times are plotted as a function of input voltage. Figure 7-1 shows 
the characteristic of TCHARG, the on time of the switching translators, 
and Figure 7-2 shows the characteristic of TDUMP, the dump time for the 
charging inductors. The off time, TDtfELL, is shown in Figure 7-3. The 
displacement of the 2 and 4 ohm load curves is due to the large currents 
flowing through the charging inductors during TCHARG. The large values 
of current causes the inductance value to decrease (approach saturation) 
resulting in a further increase in current. The outcome is current 
values which are larger than would be if the inductor did not approach 
saturation, and this in turn produces sufficient output voltage to 
require the regulator to remain in the off state longer. 

Figure 7-4 shows how the duty cycle, D, of the switching 
transistors is a function of the input voltage and the load resistance. 
Closely related to this is the information in Figure 7-5. This is the 
switching repetition rate of the regulator as a function of load resis- 
tance for two different input voltages. The shape of the curves is due 
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8.00 



V IN < volts > 


{> 10 * ) 


Figure 7-4. .Duty cycle of switching transistors as a function 
of input voltage and load res is t ance * 






LOAD RESISTANCE (ohms) 


Figure 7-5. Repetition rate of switching transistors as a 
function of load resistance and input voltage. 
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to the nonlinear decrease in TCHARG, TDUMP and TDWELL as the load resist- 
ance increases. For values of above 8 or 10 ohms, the decrease be- 
comes small for TCHARG and TDUMP while TDWELL begins to increase. 

The transient response of the regulator, with zero initial con- 
ditions, is exhibited in Figures 7-6 and 7-7. In Figure 7-6, each pulse 
of current shown is due to the charging inductors dumping their energy 
into the output filter. It is recalled that each regulator has four 
such inductors in parallel and the dump current, Iqump* will be the sum 
of their currents. The damped oscillation in I 2 , the current in the 
middle loop of the filter, is due to the natural response of the filter 
and has a frequency of approximately 250 Hz. Since the magnitude of this 
current is used to determine TCHARG, the oscillation is passed to IpxjMP 
after a small time delay. Initially, the output voltage is low enough 
to require the transistor switches to turn on immediately after the 
charging inductors have dumped. Thus the dwell time, TDWELL, is zero. 

The change in disposition of the current pulses near 10 milliseconds is 
due to the voltages of the filter becoming large enough to allow the 
regulator to dwell before requiring the switches to turn on. 

Figures 7-8 and 7-9 contain the steady-state response of the reg- 
ulator to a 4 U load and a 60 volts input. The scales of the plots have 
been expanded beyond those of the transient response to disclose more 
detail. The small change in rate of decrease of the dump current » *DUMP> 
(gee last current pulse) is due to the change in inductance of the charg- 
ing inductor. Later, another figure will examine the effect of this 
inductor on current and time values more closely. The ripple voltage 
at the input of the output filter is seen in Figure 7-9. The near 
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Figure 7-6. Transient response of filter currents. 
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Figure 7-7, 


Transient response of filter voltages. 
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Figure 7-8. Steady-state response of filter currents. 
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instantaneous jump in takes place just as the transistor switches 
begin to turn off and the charging inductors begin to dump. Its magni- 
tude is composed of two voltage terms. First , the voltage across R^, 
the effective resistance of the branch containing Cl, and, second, the 
voltage across L C1’ the effective inductance of the same branch. The 
spike component at the jump is the voltage across and the rest of 
the jump is due to the voltage across R^. As explained in Chapter III, 
the duration of the spike is equal to the turn-off time of the switching 
transistors, which is less than one microsecond. 

The end of TDUMP and beginning of TDWELL is indicated by point "a" 
on the plot, and the end of TDWELL and beginning of TCHARG is indicated 
by point 'V. The fact that the voltage across Cl doesn't appear exponen- 
tial during the discharge time, TDttEId. plus TCHARG, is due to the large 
time constant present. 

A plot of output voltage versus output current is given in 
Figure 7-10. The slope is due to the effective internal dc resistance 
of 213 milliohms programmed into the regulator by the use of current 
feedback. 

To show the nonlinear charging and discharging currents. Figure 
7-11 is now presented. This represents the worse case in that the 
regulator is working under the maximum current demand allowed, 13.5 
amperes . The value of IpuMp is the sum of the four currents coming 
from the four charging inductors. The dashed line indicates the charging 
current possible during TCHARG for the ideal case in which the inductor's 


value does not decrease. 
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The last four figures contain the response of the regulator to 
a sudden change in load resistance. Figures 7-12 and 7-13 show the 
results of the regulator, under steady-state conditions, having the 
load resistance switched from 4 ohms to 2 ohms. Figures 7-14 and 7-15 
are for the load being switched from 4 ohms to 8 ohms. 

The intent of the preceding fifteen figures was to display a 
sample of what can be accomplished with the computer program and at the 
same time familiarize the reader with some of the CSKM's characteristics. 

This report is concluded by summarizing what has been accomplish- 
ed. The necessary mathematical expressions were developed in Chapter III 
to represent the three switching times: TCHARGE, TDUMP, and TDWELL. 

Using the results of Chapter 111, a computer program model was developed 
and presented in Chapter IV. The guidelines stated in the Introduction 
were followed closely and the result was a very versatile program that 
can easily be modified to examine any particular area of interest. For 
example , in addition to what was done at the beginning of this chapter, 
it is possible to investigate the effect of varying a component value 
over a large range. 

There are at least two areas that could be investigated further. 
First , in modeling the spike present in several simplifying assump- 
tions were made. To be more exact, information is needed on the effec- 
tive inductance at the input of the output filter and a good equivalent 
model of the switching transistor is necessary. A second area is the 
calculation of the ac output impedance for very low frequencies. Due 
to the presence of nonlinearities , linear or piecewise-linear methods 
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Figure 7-14. Transient response of filter currents to a 4 ohms-to-8 
ohms load change. 
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Transient response of filter voltages to a 4 ohms-to-8 
ohms load change. 


Figure 7-15. 
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could not be used. 

Two areas that were not covered in this report are, 1) calcula- 
tions of efficiency of power transfer, and 2) an investigation of a 
damped high frequency oscillation that last for a short duration when 
the charging inductors dump. To calculate efficiency, the power loss 
in the switching transistors oust be calculated and given sufficient 
information, this can be done. To analyze the oscillation, the stray 
capacitance and inductance must be known. 
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APPENDIX A 


The development of Equation 3-2 of Section A, Chapter III, re- 
peated below as Equation A-l, will be done by solving for 1^. 


! a (S) = 


-V, 


a 3 S 3 + a 2 S 2 + a^S 1 + a Q 


V \l + + *L S b 4 S 4 + b 3 S 3 + b 2 S 2 + ^S 1 + b Q 


(A-l) 


The form of I and I _ will be as that of i.,. The matrix expression 
A2 A3 A1 

for the filter currents was given in Equation 3-1 and using Cramer's 
rule for the solution of equations will yield: 


a 4 s 4 + * • • + a 1 s-v d 

S(b 4 S 4 + - • + bjS + b Q ) 

The above equation can be expanded into two terms. 

Iil<S) - ^ + — 7 

s b.S 4 + **• + b.S + b rt 
4 10 

Solving for and K 2 yields: 


I A1< S > = 

** D 


(A-2) 


(A-3) 
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*1 



(A-4) 


and 


K2 " (A^S 3 + • * • + A^) + (b 4 S 3 + • * * + b ± ) (A— 5) 

It Is noticed that the first term of K2 Is the numerator of D^/D 
with the last term deleted and the powers of S decreased by one. In 
addition, the second term of K2 is V^/b^ times the denominator of D^/D 
with the last term deleted and the powers of S decreased by two. Using 
this approach, the K^'s are calculated and called CNUMA in the computer 


program. 



APPENDIX B 


In this section, the parameters of Equation 3-4, repeated below, 
are discussed . 

ia(t) “ 2R a1 expt-a^t] cos(wyt + 8^) 

+ expt-a^t] cosCto^t + 9^) 

“ V d /(R d + R L1 + R L 2 + Rj,) (B-l) 

This equation represents the form of the three loop currents of the out- 
put filter while it is in mode A. The notation used below is the same 
as that used in the computer program. 

The three currents are called CURNT(l), CURNT(2) , and CURNT(3) 
in the computer program and each subscript, (1), (2), and (3) refers to 
a particular loop in the filter. The method of the Heaviside partial 
fraction expansion can be programmed to determine the parameters of these 
three currents. First, there are six RA's to be found which are the 
magnitudes of the coefficients of the partial fraction expansion. They 
will be double subscripted showing first the current to which they belong 
and second showing position within that current equation. Thus the six 
RA's are: 
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RA(1,1) > RA(1,2) 

RA(2,1) , RA(2 ,2) 

RA(3,1) , RA(3,2) 

In general, the magnitude RA^ for a complex root can be found from: 

CKA = RA Ak ± - (S + oA* - juAj) Jfff < B “ 2 ) 

1 - L Q(S)j s * -ctAi + jwA* 

Using this approach for example, R(l,l) will be: 

|CKA(1,1)| - RA(1,1) - |- C S^| (B-3) 

where CNUMA(l,k) is the numerator of the polynomial fraction in I^(S) 
evaluated at S equal to -aAl + jwAl. DNOMA(k) is the denominator of the 
polynomial fraction in I^(S) less the S + aAl - joiAl factor, evaluated 
at S equal to -aAl + jaiAl. 

The ciA's, called ALPHAA(l) and ALPHAA(2) in the program, are 
equal to minus the real parts of the complex roots used to determine the 
RA's. The wA's, called CMEGAA(l) and 0MEGAA(2) in the program, are the 
positive imaginary parts of the complex roots used to determine the RA's. 
Note that each current equation will contain the two aA's and toA's. 

The 6A*s are the phase angles associated with the coefficients 
of the partial fraction expansion that were used in determining the RA's. 
They will be double subscripted showing first the current to which they 
belong and second showing position within that current equation. 
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The six 0A's are 

THETAA(1,1) t THETAA(1,2) 

THETAA(2 , 1) , THETAA(2 ,2) 

THETAA(3,1) , THETAA(3,2) 

They can be found by taking the arctangent of the imaginary part of the 
coefficients used in determining the RA*s divided by its real part. 

For example : 


THETAA(1,1) - Arctan 



(B-4) 


The procedure used in finding the parameters for the mode A 
current equations is used to find those belonging to the equations for 


mode B. 



